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Abstract 

In this paper we apply variational energy band theory to a form of the Holstein 
Hamiltonian in which the influence of lattice vibrations (optical phonons) on both local 
site energies (local coupling) and transfers of electronic excitations between neighboring 
sites (nonlocal coupling) is taken into account. A flexible spanning set of orthonormal 
eigenfunctions of the joint exciton-phonon crystal momentum is used to arrive at a 
variational estimate (bound) of the ground state energy for every value of the joint 
crystal momentum, yielding a variational estimate of the lowest polaron energy band 
across the entire Brillouin zone, as well as the complete set of polaron Bloch func- 
tions associated with this band. The variation is implemented numerically, avoiding 
restrictive assumptions that have limited the scope of previous assaults on the same 
and similar problems. Polaron energy bands and the structure of the associated Bloch 
states are studied at general points in the three-dimensional parameter space of the 
model Hamiltonian (electronic tunneling, local coupling, nonlocal coupling), though 
our principal emphasis lay in under-studied area of nonlocal coupling and its interplay 
with electronic tunneling; a phase diagram summarizing the latter is presented. The 
common notion of a "self-trapping transition" is addressed and generalized. 
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1 Introduction 

In a previous paper [p], we undertook a numerical implementation of methods proposed by 
Munn and Silbey || for the determination of polaron properties in the presence of simultane- 
ous local and nonlocal exciton-phonon coupling, local coupling being defined as a nontrivial 
dependence of the exciton site energies on lattice coordinates and nonlocal coupling as a 
nontrivial dependence of the exciton transfer integral on lattice coordinates. Those meth- 
ods were essentially perturbative, with the coefficients of a canonical transformation being 
fixed so as to limit the secular growth with the temperature of the perturbation remain- 
ing after transformation. The present work is motivated by several limitations of those 
methods that were found in the course of our prior work [|TJ: 1) Since the Munn-Silbey 
method was motivated by an interest in controlling perturbations at high temperatures, one 
should not necessarily expect the method to produce the best description of polaron states 
at low temperatures. 2) Since the method is perturbative, it is reasonable to expect that 
the description of the polaron ground state so obtained could be improved. 3) Since the 
self-consistency equations on which the method is based do not involve the exciton transfer 
integral, it is questionable whether the method can be relied upon beyond the very narrow 
band regime. 4) Since the results of our numerical investigation showed some significant 
differences with analytical results based on the same underlying methods, as well as with 
some other recent analyses 0-0, it would appear desirable to have corroboration of our 
results by independent methods. 

The central findings of our prior work [ij, that polaron binding energies may be much 
larger than previously thought, and polaron bands much narrower, underscores the impor- 
tance of independent corroboration, since such findings are of central importance to under- 
standing the influence of nonlocal exciton-phonon coupling on the nature of polaron states 
and of polaron transport. 

Simultaneous local and nonlocal coupling appears to be particularly important in the 
characterization of solid-state excimers, where a variety of experimental and theoretical con- 
siderations suggest that a strong dependence of electronic tunneling upon certain coordinated 
distortions of neighboring molecules (nonlocal exction-phonon coupling) is crucial to the for- 
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mation of excited bound states g § 0, §, g 0, 0, [12|, |I§ 0, |I§ [Tj| |T7|, |T|, [19], [§, [21], [22 



23, 24, 25 



In this paper we approach the problem of simultaneous local and nonlocal exciton-phonon 
coupling by variational methods. 

Our central interest in this paper is in the polaron energy band, computed as 

E K = (V!(k)\H\V(k)) (1) 

wherein k is the total crystal momentum label, ^(ac)) is an appropriately normalized delo- 
calized trial state, and H is the system Hamiltonian. It should be noted that the trial states 
we use are eigenf unctions of the appropriate total momentum operator, so that variations for 
distinct k are independent, and the set of E K so produced constitute a variational estimate 
(upper bound) for the polaron energy band |26|, [27]. 

As our system Hamiltonian, we choose perhaps the simplest one embracing exciton tun- 
neling and simultaneous local and nonlocal exciton-phonon coupling, a slight generalzation 



of the traditional Holstein Hamiltonian [28 



H = H ex + H ph + H ex - ph (2) 

^ ex, = -«/E4K+i + «n-i) (3) 

n 

W h = coJ2b% (4) 

n 

fyex-ph _ gu^alfinfy + b n ) + -U Y^{4i a n-1 + i-l^n ~ ^+1°" _ 4 a n+l)(^ + K) (5) 



in which cv n creates an exciton in the rigid-lattice Wannier state at site n, and b* n creates 
a quantum of vibrational energy in the Einstein oscillator at site n. Both excitons and 
phonons are treated as bosons. The Einstein frequency is given by u, J is the exciton transfer 
integral between nearest neighbor sites, g is the local coupling strength, and </> is the nonlocal 
coupling strength characterizing phonon assisted transfers between nearest neighbor sites. 
All of the methods used in this paper can be applied as well to common generalizations 
of this Hamiltonian, involving, for example, phonon dispersion, exciton transfers beyond 
nearest neighbors, and/or different exciton-phonon coupling geometry. 
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In Ref U] we considered both antisymmetric nonlocal coupling, as above, and symmet- 
ric nonlocal coupling, in which all nonlocal coupling terms would be of the same algebraic 
sign. These two kinds of coupling represent different physical circumstances. Antisymmetric 
coupling, for example, would be appropriate to the description of certain librations that 
promote exciton transfers between neighboring molecules when these molecules tilt toward 
each other, effectively decreasing the gap through which tunneling must take place. This 
has the consequence that tunneling between a molecule and its neighbor to the right (for 
example) is promoted (and tunneling on the left inhibited) when the librator tilts to the 
right, and tunneling on the left is promoted (and tunneling on the right inhibited) when the 
librator tilts to the left. Symmetric coupling, on the other hand, describes the circumstance 
in which tunneling between a molecule and its neighbors on both the left and right is pro- 
moted during the same phase of oscillation and inhibited during the complementary phase; 
this may happen, for example, if the strength with which a mobile exciton is bound varies 
with the coordinate described by the oscillator. In this paper, we restrict our attention to 
antisymmetric nonlocal coupling. 

In momentum space, these several terms take the form 

H ex = J2 J k a\a k (6) 

k 

H* h = Y,ub\b q (7) 
i 

H ex - ph = iv- 1/2 E "/ViW&* + fet -«) ( g ) 

kq 

where 

Jk = —2J cos k, fl = g + z^[sin k — sin(k — q)], (9) 

Throughout this paper, we use the the Fourier conventions for ladder operators c' = a\b\ 
etc. and scalars 7 = a, {3, etc. 

4 = iv" 1/2 E z~ ip % 4 = N ~ 1/2 E ( 10 ) 

p n 

ln = A^E^Tp, 7 P = E^ pn 7n, (11) 

p n 
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To assist in the interpretation of the various formulas to follow, we denote exciton wave 
vectors by latin Fs, phonon wave vectors by latin g's, and reserve the greek n for the total 
crystal momentum label. 

Also throughout this paper, we use the Einstein frequency as the scale of the energy, so 
that all energies are dimensionless. This leaves us with a three-dimensional parameter space 
in which a general point may be identified as (J, g, <fi). 

A complete exploration of the (J,g,(p) parameter space is far beyond the scope of this 
paper, so we must narrow our focus significantly. Since our principal interest here is in 
nonlocal exciton-phonon coupling, we defer an in-depth discussion of the (J,g,0) plane to 
other work. Since we wish to compare our present results with those of the Munn-Silbey 
method, we present results for the (0,g,(p) plane in a manner that closely follows prior 
work. Since the regime of strong nonlocal coupling has received much less attention than 
that of local coupling, we present a more systematic survey of the (J, 0, <fi) plane, including a 
phase diagram and a discussion of characteristic polaron structures and the nonlocal-coupling 
version of the self-trapping transition. Finally, we present sample results for general points 
in the (J, g, (/)) space. 

The layout of the paper is as follows. As a preliminary study, we briefly touch upon 
the application of the small polaron Ansatz to nonlocal exciton-phonon coupling, as was 
originally suggested by Merrifield in his 1964 paper |5C|| , but was never carried out. Toy- 
ozawa's Ansatz is then introduced to examine the same problem in greater depth. The rest 
of the paper is devoted to exploring the parameter space under Toyozawa's Ansatz, and 
a comparison of our present results with comparable results of the Munn-Silbey approach 
is made. We must emphasize that these comparisons are not to the original analytical 
calculations of Munn and Silbey 0, but to our own implementation of their approach by 
numerical methods . Since the results of our numerical study differed quantitatively and 
qualitatively from the original calculations of Munn and Silbey in some significant respects 
and thereby generalized them, it is important in the following to distinguish between the 
Munn-Silbey method and our own numerical results obtained by this method. 
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2 The Small Polaron Ansatz 

As in Merrifield's original calculation, the normalized small polaron trial state may be written 

!*(«)> = i^ 1/2 E eiKn 4ex P [-E(^ 2 -n^ 2 -C-AJ]|o>, (12) 

n n2 

|*(«)> = iV-V^ e i- a t exphiV -i/2^ (/3 « e ^n & t _ ^V«%)]|0). (13) 

n g 

This trial state may be viewed as a phased sum over "form factors" describing local exciton- 
phonon correlations. The form factor in this case is the product of an exciton completely 
localized on a single lattice site and a lattice function defined relative to the exciton. 

Using this trial state, we form the expectation value of the Hamiltonian and minimize 
the total energy E K with respect to the phonon amplitudes, yielding the self-consistency 
equations 

= g + ^^cos^-^-fjsinf 

Pq 1 - 4JS"mn(K -$*-§) sin f - 805^ sin § ' 1 ' 

5* = aq)[JV- 1 X;i^| a (ooBg-l)] J (15) 
$ K = AT- 1 Easing, (16) 

9 

^ = AT 1 E Mfl?) ™(« - $ K - | - |) sin Q -. (17) 

9' 

in which S K and $ K are the magnitude and phase of the Debye- Waller factor (see Section 
III). 



It is evident from ( fn| ) that unlike the original work of Merrifield restricted to local 
coupling only f30fl , (3g here can not be taken to be real owing to the presence of nonlocal 
exciton-phonon coupling. The phonon mode amplitude f3* is real for all k and q when = 0; 
when (ft 0, (3* is real only along the lines defined by [sin(fc — $ K ) — sin(« — $ K — q)] = 0, 
among which is the q = line. 

Using a numerical iteration scheme, we may solve the above set of self-consistency equa- 
tions to desired precision without any artificial constraints. For all values of local and non- 
local coupling strengths, we found our final solutions to be insensitive to our choice of both 
numerical method and states used to initialize calculation. For the purpose of illustration, 
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we may express (3* in terms of two real matrices 77* as in Ref 

0," = + ^<[sin(K - $ K ) - sin( K - $ K - g)]. (18) 

From Eq. (|TJ|) it follows that 77" = S K . The (/-independence of S K suggests that the £ and 
77 matrices may be more similar in shape than we found in our numerical implementation of 
the Munn-Silbey method |lj; moreover, since K-dependence of S K can be relatively weak in 
some regimes, the difference between £ and 77 can be quite small. 

Results of such calculations for weak nonlocal coupling (0, 1, 0.1) are indicated in Fig. |l[ 
Only £^ is displayed since S K in this case varies by only a fraction of a percent across the entire 
Brillouin zone so that 77^ would be indistinguishable from Since nonlocal coupling is weak 
relative to local coupling in this particular example, it follows as well that (3^ ~ The 
weak asymmetry of (3* with respect to q at intermediate values of k implies (through 
Eq. [IB]) that the Debye- Waller phase $ K varies slightly with k; however, this weak non- 
constancy of $ K has little direct influence over the properties of the solution both because 
of the smallness of $ K and the small value of the nonlocal coupling constant. 

Antisymmetric nonlocal exciton-phonon coupling introduces characteristic distortions in 
quantities such as and the polaron energy band. These distortions are "bimodal" with 
respect to the total crystal momentum k in the sense that this modulation characteristically 
is most pronounced at intermediate values of \k\. With respect to the phonon wave vector, 
these distortions characteristically are most pronounced at |g| = n, indicating that the 
strongest exciton-phonon coupling involves lattice dimerization. On the other hand, ^ is 
weakly structured at low phonon wave vectors, and is in fact equal to unity along q = 0, 
consistent with the lattice "sum rule" 

(3 q=0 = J2Pn = g ■ (19) 

n 

The polaron band energy E K can be calculated from 

E K = N-^W -2JS K cos{k-<5> K )-2N- 1 gY,Re{Pq ) 
g 1 

-an- 1 4>s* Im (^) cos (^ - $K - h Mb- (20) 

q 11 

Comparison between the small polaron band and that calculated numerically from the Munn- 
Silbey approach [|IJ is shown in Fig. |2] for the case of (J,g,(j)) = (0, 1, 1). 
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Figure 1: The phonon displacement factor calculated from the small polaron Ansatz for 
(J, g,4>) = (0, 1, 0.1). The variation in is small in this case because the nonlocal coupling 
strength </> is weak relative to the local coupling strength g. 

The Munn-Silbey method yields a flatter and lower energy band than does the small- 
polaron method; however, since the perturbative result falls below the variational result, no 
conclusion can be drawn as to which approach gives the better polaron band; for example, 
without further information we cannot exclude the possibility that the actual ground state 
band might lie between these two results and closer to the variational band. Higher-quality 
variations are needed to make such determinations. 

Through its Fourier transform, the nonuniformity of with respect to q (though small 
on an absolute scale in this case) implies a spreading of phonon distortion in site space 
even in the absence of a direct transfer integral J. With increasing nonlocal coupling, the 
amplitude of this distortion increases, as does its spread in real space. 

Since nonlocal exciton-phonon coupling provides the only transport mechanism when 
J = 0, it is reasonable to expect in analogy with the usual small polaron picture (e.g., as 
in Ref. fl30| ) that a small polaron approach such as this one may break down when nonlocal 
coupling becomes comparable in strength to local coupling. 
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Figure 2: Comparison of the polaron bands between the Munn-Silbey approach |l[ (points) 
and the small polaron Ansatz (solid line) for (J,g,(p) = (0,1,1,) at T = 0. The small 
polaron Ansatz gives a higher ground state energy and a larger bimodal variation in the 
polaron band. 

At the least, the approach sketched here needs to be generalized to admit wave function 
symmetries more consistent with the dimeric structures favored by nonlocal coupling. This 
could be done in a minimal way by building up trial Bloch states from dimeric form factors; 
we elect, however, to skip this step and proceed to a still more flexible trial state introduced 
by Toyozawa that is well known in the context of local-coupling polaron theory. 

3 Toyozawa' s Ansatz 

The major shortcoming of the small polaron Ansatz is the fact that it is built up from form 
factors in which the electronic component is completely localized on a single lattice site. 
In the presence of significant exciton tunneling, or, as we have observed, in the presence 
of significant nonlocal exciton-phonon coupling, it would appear reasonable to expect self- 
consistent local exciton-phonon correlations to reflect at least some of the characteristics of 
spreading wave packets. To accommodate this expecatation, the completely-localized exciton 
component of the small polaron form factor can be generalized to a superposition of exciton 
amplitudes to be determined self-consistently with the amplitudes of the lattice oscillators. 
This is the basis of many theories characterizable as "large polaron" theories. Important 
among these for our purposes are delocalized-state theories based on phased sums of such 
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generalized form factors (historically identified with Toyozawa), and localized-state theories 
based upon using such generalized form factors as trial states in their own right without 
invoking the phased-sum construction that assures derealization (historically identified with 
Pekar 0§|,§ and Davydov |5|, |3|, |37|, g gg, gl], 0, g|, g|, gl| gl], g7[). 



Toyozawa's Ansatz may be viewed as a derealization of the Davydov Ansatz |?/>) 

\1>) = |a>®l£>, (21) 

where Cg> denotes the direct product, and \a) and \(3) are the exciton and phonon part of the 
form factor respectively, built from exciton and phonon vacua \0) ex and (0)^ as 

|a) = X;on4|0)ex, (22) 

n 

\P) = ew\EM-ftK)]\0)»- (23) 

n 

After derealization one obtains Toyozawa's Ansatz state, given by 

= (24) 

|«> = E eJKn <-n<exp[-E(^ 2 -n^ 2 -PT 2 -nKW) (25) 

nni ri2 

|k) = N- 1/2 J2e i{K ~ k)n a K k alexp[-N- 1/2 Y,((3ge- iqn bl - p«*e iqn b q )]\0). (26) 

nfc g 

where |0) = |0} ex (8) \0) p h- The auxiliary vector |k) is not normalized, but simplifies presenta- 
tion of some results. It may be worth noting that many localized-state approaches such as 
the Davydov theory take the exciton component |a) to be normalized such that J2n \ a n\ 2 — 1) 
since in this approach the norm of exciton amplitudes has a physical interpetation as the 
exciton number. In the delocalized-state approach, however, the norm of the exciton ampli- 
tudes does not have this physical meaning, and by contributing only a multiplicative factor 
to makes no contribution to the normalized trial state. Consequently, the norm of the 

amplitudes a is arbitrary, and in fact drifts in the course of numerical variation; since this 
drift has no physical consequences, we have found it convenient throughout this paper to 
normalize the exciton amplitudes such that a n= o = 1. 

For variational calculations we require the total energy associated with each trial state, 
for which we evaluate the expectation values of the three principal terms in the Holstein 
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Hamiltonian: 

(K\H ex \K) = -2JN- 1 S%-k cos k K\ 2 , (27) 

k 

(k\H»\k) = iV- 2 E^_ fc _ 9 |« fe 1 2 |/3 g 1 2 , (28) 



kq 



h^\k) = -N- 2 y: a- q <at +q (s:_ k _ g %* + s:_ k ^ q ). (29) 

kq 

(^) = N- 1 J2s:-kK\ 2 - (30) 

k 

Here, S% is the Fourier transform of the generalized Debye- Waller factors S% 

^ = Ee-^, (31) 

n 

SZ = e X p[N- ly E\j3*\*(t* n -l)], (32) 

which quantifies the overlap between the lattice components of polaron wavefunctions dis- 
placed from each other by n lattice sites. These are to be distinguished from the Franck- 
Condon factor 

ph W) = e~^n\M 2 , (33) 

which quantifies the overlap between the lattice component of the polaron wavefunction with 
the undistorted phonon ground state. The Debye Waller factors <S^ =±1 (= S K e ± ' t ' s ' K of section 
II), appear routinely in the transport terms of effective (small) polaron Hamiltonians, where 
they strongly influence the renormalization of the effective mass. In the general case we 
address here, the spread of the exciton amplitudes a n causes Debye- Waller factors between 
non-nearest-neighbor sites to contribute to polaron structure in complex ways. 

Owing to the large number of Debye- Waller factors that may be involved in the general 
case, their resolution into magnitudes and phases as in section II ceases to be advantageous. 
One can show from ( |3"T| ) and ( |3"2] ) that S k is strictly real, so that computation is simplified 
by abandoning the magnitude and phase characterization and working in momentum space 
with the real quantities S%- 

Minimization of the total energy E K with respect to the phonon amplitudes (3** leads to 
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where 

L K q = N- 1 £ a_ q S:_ k _ q aTa K k+q , (35) 
M^iV-^S^J^I 2 , (36) 



and H q is the sum of three terms: 



H? = -2JN- 1 '£Sl. k _ q <x*k\aZ\*, (38) 
k 

Hf = N- 2 J2S K K _ k _ q _ q ,\at\ 2 \P«,\ 2 , (39) 

kq' 

H e q x ~ ph = ~ N ' 2 £ + ^*_ 9 ^). (40) 

Similarly, after minimizing with respect to a k * we arrive at 

= Afjf — (E" + 2Jcosjfc)^_ fc ^ 

where 

= N- 1 £ fZL q al +q {S^ k _ q (i^ + S:_ k (31 q ) (42) 

and 

M^iV- 1 ^^^! 2 . (43) 
i 

We may recover the lattice sum rule (19) from (|34]) by noting that 



H* =0 = £?«(«|«> , ^L^,, = M g K =0 = («|«>. (44) 

4 p and 



The effect of nonlocal exciton-phonon coupling is most evident when J is small or absent. 
Moreover, since local and nonlocal coupling influence polaron structure in distinct ways, it 
is well to examine the influence of nonlocal coupling both in isolation and in concert with 
local coupling. 

The first panel of Fig. ^| shows the exciton and phonon amplitudes (3% and a 1 ^ for the 
pure nonlocal coupling case (J,g,(f>) = (0,0,1) at k — where both of these quantities are 
real. 
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Figure 3: Variational parameters for the k = state calculated from Toyozawa's Ansatz. 
Second panel: the exciton amplitude a* =0 (dashed line) and the phonon displacement /3^ =0 
(solid line) for (J, g, (/)) = (0, 0, 1); First panel: the exciton amplitude a* =0 (dashed line) and 
the phonon displacement /3^ =0 (solid line) for (J,g.<p) = (0,0.3,1). 

In the absence of local coupling, the exciton amplitudes o^ =0 are bond-centered and 
even under inversion through the occupied bond, while j3* =0 is bond-centered and odd under 
inversion, reflecting the antisymmetry of the nonlocal exciton-phonon coupling. Approx- 
imately 84% of the total exciton density is shared equally by the two sites defining the 
central bond. This is characteristic of the strong nonlocal coupling regime, and justifies 
some the assumptions of our previous work. 

The first panel of Fig. should be compared to the exciton-phonon correlation function 
Ci produced as a diagnostic of polaron structure in our previous paper |1J; C\ was constructed 
around a dimeric exciton function ip n restricted to two sites (i.e., 100% of the total exciton 
density shared equally between n = 0, 1) that we presumed to be representative of the exciton 
amplitudes a* =0 in actual polaron states. The similarity between (3^ =0 of Fig. ^]a and Ci is 
striking. (The previous calculation included a small local coupling strength that does not 
significantly affect the comparison.) 

The second panel of Fig. |^ shows (3^ =0 and a* =0 for (J,g,<p) = (0,0.3,1.0), reflecting 
the addition of a moderate amount of local coupling to the previous scenario. It is still the 
case that most of the total exciton density resides on the two sites defining the central bond 
(this does not change materially until J becomes signficant); however, the addition of a non- 
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negligible amount of local coupling destroys the bond-centered symmetry by introducing a 
site-centered structure localized on only one of the two central sites. 

Fig. f| displays the energy band comparison between the Munn-Silbey approach |l| and 
Toyozawa's Ansatz for (0,0.03,1) (essentially Figure |a) and (0,1,1), both of which were 
studied in Ref. 0. 
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Figure 4: Comparison between the polaron band calculated from the Munn-Silbey approach 
Hat T = 0, (diamonds, (0, 1, 1), crosses, (0, 0.03, 1) and that from Toyozawa's Ansatz (solid 
line, (0, 1, 1), dashed line, (0, 0.03, 1). 

Though still imperfect, the energy band comparison is much more favorable here than we 
found in section II. The consonance between our present results and our previous calculations 
by the Munn-Silbey method appears to be due to the ability of both approaches to embrace 
dimeric exciton-phonon correlations essential to nonlocal-coupling polarons. Our present 
results based on Toyozawa's Ansatz give the best estimate of the polaron ground state 
energy among the three we consider, and can be viewed as corroborating our prior numerical 
calculations by the Munn-Silbey method, at least when nonlocal coupling dominates. Due 
to the variational nature of Toyozawa's Ansatz, it is now safe to say that the Munn-Silbey 
approach over estimates the ground state energy at T = 0. 
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5 J and <\> 

In this section we put local coupling aside (i.e., we set g = 0), and focus on the interplay 
between exciton tunneling (J) and nonlocal coupling (<fi). We have seen that although 
nonlocal coupling is an exciton-phonon interaction that supports polaron formation, it is 
also a transport mechanism and competes with local coupling both by promoting transport 
and by driving the exciton-lattice correlations toward dimeric structures rather than the site- 
localized structures preferred by local coupling interactions. On the other hand, although 
nonlocal coupling is a transport mechanism, it also competes with direct, phonon-free exciton 
transfers since the lattice distortions inherent in phonon-assisted transfers inhibit direct 
transfers. 

It is convenient to use the device of a phase diagram to organize our discussion. Fig |5] 
depicts the J — plane at g = 0, in which we have indicated two narrow, tongue-shaped 
regions that divide the plane into two more-or-less distinct regions, a strong nonlocal coupling 
region occupying the upper portion of the diagram and a weak nonlocal coupling region 
occupying the lower. The phenomena indicated by these tongues can be characterized as a 
kind of self-trapping. 

6 
5 
4 

Q. J 

2 
1 


2 4 6 8 10 

J 

Figure 5: Phase diagram on the J — <fi plane for Toyozawa's Ansatz. The two wedges 
correspond to the discontinuity near the Brillouin zone center (upper wedge) and that near 
the Brillouin zone boundary (lower wedge), respectively. 

In the more familiar self-trapping phenomenon associated with local exciton-phonon cou- 
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pling, a polaron is understood to be characterized by exciton-phonon correlations that are 
more-or-less spatially compact. For "weak" local coupling, correlations may extend over 
many lattice sites, corresponding to a "large" polaron. With increasing local coupling 
strength, the width of the correlated region diminishes until the excitation is essentially 
confined to a single lattice site, corresponding to a "small" polaron. Along the way, other 
polaron properties change as well, most notably the polaron effective mass growing from 
the free mass at weak coupling to values arbitrarily large in the strong coupling region. 
Whether this change as a function of system parameters is necessarily smooth or may occur 



discontinuously is a subject of some contention [48, 49, BQ, 51]; nonetheless, it is common 



for approximate treatments to find the change to be discontinuous in some regimes, and 
hence the imprecise notion of a discrete "self-trapping transition" is widespread. More im- 
portant than the question of continuity is the fact that in the vicinity of this transition, 
exciton-phonon correlations and related physical quantities such as the effective mass un- 
dergo strong changes with relatively small changes in control parameters. 

The qualitative characteristics of this familiar self-trapping transition apply as well to 
our treatment of nonlocal coupling. We find a discrete transition to occur in the vicinity 
of the Brillouin zone center, bearing some resemblance to the more familiar local coupling 
phenomenon; however, we also find a discrete transition occurring in the vicinity of the zone 
edge. Moreover, whereas the usual conception of self-trapping focusses on dramatic changes 
that occur at the Brillouin zone center (e.g., the jump in the effective mass) leading to the 
notion of a sharp transition, we find the discrete transitions in both the inner and outer 
zones to be broad in a sense to be described presently. 

We must now attempt to be precise in specifying the meaning of "discrete transition" in 
the present context, for clarity focussing on the inner-zone transition. 

At an arbitrarily chosen point (J, g, 4>) in parameter space, there exist one or more minima 
in the variational energy E K for each k, the lowest of which identifies the global energy 
minimum in each total momentum sector; the collection of energies for all k identifies the 
polaron energy band E K at (J, g, </>), and the collection of states associated with these global 
energy minima constitute the polaron Bloch states. The distinct classes of relative minima 
that may coexist over certain invervals of k are associated with distinct classes of variational 
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states corresponding to different polaron structures. In the present problem, typically, either: 
1) one class of minima exists for all k, or 2) one class of minima exists over an interval 
\k\ G (0, Ki), and another exists over an interval \k\ G («2, tt), with the two classes of relative 
minima coexisting for \k\ G (k 2 , ki). In the latter case, there exists a particular total crystal 
momentum n* G («2,Ki) at which the polaron energy band changes from being defined by 
one class of global energy minima in the inner zone < k*) to being defined by the other 
class of global energy minima in the outer zone > k*); this k*, like the energy band 
itself, is a function of the system parameters ( J, g, (f>). 

We describe the 1ow-ac states as "large-polaron-like" and the high-K states as "small- 
polaron-like" because through sequences of infinitesimal steps in k and/or exciton-phonon 
coupling parameters it is possible to smoothly deform any large-polaron-like state into a 
traditional large polaron state at n — 0, and any small-polaron-like state into a traditional 
small polaron state at k — 0. Significant for the explication of the discrete transition phe- 
nomenon is the fact that it is not generally possible to smoothly deform large-polaron-like 
states into small-polaron-like states (or vice versa) through such sequences. In particular, at 
fixed (J,g, 4>) smooth changes in k are followed by smooth changes in all polaron properties 
until k reaches k*. At k*, discontinuities appear in at least some polaron properities, conse- 
quent to the switching of the global energy minimum from one class of state in one part of 
the variational state space to another class of state finitely separated from the first. 

The appearance of a ^-dependent discrete transition at n* is unambiguous in our results, 
and by our description is a discrete transition between small- and large-polaron-like states. 
What may be less clear to this point is what relationship, if any, exists between this transi- 
tion phenomenon and the more traditional notion of a discrete self-trapping transition. The 
"order parameter" commonly used as an indicator of the traditional self-trapping transition 
is the polaron effective mass, a quantity defined at n = 0; as such, its characteristic jump 
indicates for us the point in parameter space at which the k = polaron Bloch state switches 
from being large-polaron-like to small-polaron-like. It is possible to connect this traditional 
perspective of a n = transition with our present ^-dependent transition through the de- 
pendence of k* upon the system parameters (J,g,<f>). The general result is that whenever 
a ^-dependent transition exists between small- and large-polaron-like states at a finite k*, 
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there exist sequences of infinitesimal changes in the parameters (J,g,4>) that cause k* to 
vanish. This vanishing of k* has the consequence that the large-polaron-like region ceases 
to exist, rendering the entire polaron band small-polaron-like, including the k = state. 
This shows the usual self-trapping line to mark a boundary of the finite region of parameter 
space over which the K-dependent transition exists. This smooth connection also allows us 
to characterise the ^-dependent transition as a continuation to finite k of the traditional 
self-trapping transition at k = 0. 

Where such finite-K transitions occur, some surgery must be performed to arrive at the 
polaron band that at every k identifies the global energy minimum. This surgery consists 
of discarding the higher-lying of the coexisting variational energy minima in the interval 
(k 2 , k±) and the states associated with them, retaining only those yielding the lowest energy 
at every n. The result is a polaron band that is large-polaron-like for \k\ < k* and small- 
polaron-like for \k\ > k* . Outside of the transition regions, polaron bands are determined 
without necessity of surgery; however this qualitative characterization of the k dependence 
of the polaron band as being more large-polaron-like at long wave lengths and more small- 
polaron-like at short wave lengths continues to apply. 

The outer-zone transition is qualitatively similar in most respects, except that instead 
of going to zero in some limit, the k* associated with the transition can be pushed to the 
Brillouin zone edge, and the polaron states on each side of the transition differ from those 
of the inner- zone transition in some characteristic ways addressed later in this section. 

The tongue-shaped regions indicated on the phase diagram in Figure |^ identify param- 
eters associated with such discrete transitions; the upper, wider tongue is associated with 
the inner-zone transition, and the lower, more narrow tongue with the outer-zone transition. 
Where the tongues overlap, discrete transitions are found to occur in the inner and outer 
zones simultaneously, such that the same energy band embraces distinct polaron states in 
the inner, intermediate, and outer Brillouin zone. 

The way in which these discrete transitions depend on the parameters J and <fi can be 
sketched qualitatively as follows: For fixed J sufficiently large, a finite k* first appears as (f) is 
increased from zero to the weak-coupling boundary of the zone-center transition region. As 
is further increased, this finite k* decreases until at the strong-coupling boundary of the 
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zone-center transition k* vanishes, identifying the strong- coupling boundary with the usual 
notion of a discrete self-trapping transition occurring at the zone center. Transits across the 
outer-zone transition region have a similar structure. 

We can study the nature of the "large" and "small" polaron states associated with non- 
local coupling by examining the structure of typical states at selected points in the system 
parameter space. 

We turn first to the inner-zone transition and specifically the point (J,g,<p) = (6,0,4). 
This point is near the strong-coupling edge of the transition region, where two distinct classes 
of stable states exist in a small interval around k = 0. The structure of the two stable k = 
solutions is illustrated in Fig. ||. 



2.5 




Site Number Site Number 



Figure 6: Two types of convergent solutions for Re[[3* =0 } (first panel) and i?e[a^ =0 ] (second 
panel). (J,g,4>) = (6,0,4). The solid line is obtained when this point is approached from 
the strong coupling regime, and the dashed line is obtained when this point is approached 
from the weak coupling regime. 

The solid lines are iterative results obtained when this point is approached from the 
upper portion of the phase diagram. The dashed lines are results obtained upon approach 
from the lower portion of the phase diagram. Consequently, the solid lines show the polaron 
structure typical of the upper portion of the J — <p plane where nonlocal exciton-phonon 
coupling is relatively strong, and the dashed lines show the polaron structure typical of the 
lower portion of the phase diagram where exciton-phonon coupling is relatively weak. As 
seen in Fig. ^ the strong-coupling polaron state (solid curves) are characterized by a larger 
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phonon displacement /3^ =0 and a slightly more localized exciton amplitude a£ =0 than those 
of the weak-coupling polaron state (dashed curves), consistent with the notion of small and 
large polarons, respectively. 

The "size" of the polarons represented by the two solutions in Fig. ^ do not differ 
markedly, although upon close observation the exciton and phonon amplitudes of the "small 
polaron" state (solid curves) do decay somewhat more rapidly than do those of the large 
polaron state (dashed curves). A feature more characteristic of the difference between the 
large and small polaron states in the presence of nonlocal coupling is the fact that amplitudes 
in the small polaron state are characterized by more or stronger alternations of sign than 
are found in the large polaron state. Although these alternations are ultimately due to the 
antisymmetric nature of the nonlocal coupling used in this paper, their greater prominance 
on the strong-coupling side of the zone-center transition is consistent with the more compact 
exciton-phonon correlations typical of small polarons. 

Next we turn to the outer- zone transition, and specifically the point (J, g, 4>) = (6, 0, 3.2). 
This point is near the strong-coupling edge of the transition region, where two distinct classes 
of stable states exist in the neighborhood of \k\ = ir. Structures typical of these two kinds 
of solution are illustrated in Fig. [7[ 

As in Fig. || the solid curves in Fig. |7| are obtained when this point (J, g, (f>) = (6, 0, 3.2) 
is approached from the strong-coupling regime, while the dashed curves are obtained by 
approaching from the weak-coupling regime. A significant feature of the zone-edge solutions 
is that the phonon amplitudes can be resolved into two components; one component is a 
strongly localized asymmetric lattice distortion as we have seen in all the other nonlocal 
coupling scenarios, and the other component is a nearly-uniform plane wave. The zone-edge 
transition can be viewed as a binding or unbinding of the free-phonon component suggested 
by this plane wave. In the large polaron regime below the transition (higher k, weaker 
coupling), the plane- wave component dominates the phonon amplitudes, while after crossing 
over into the small polaron regime (lower k, stronger coupling), this plane- wave component 
is absent. 

In the particular solution illustrated, the plane-wave component has the wave vector 
q = ir, consistent with the notion of a zone-edge solution in which the phonon component 
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Figure 7: Two types of convergent solutions for Re[(3^ =n ] (first panel) and Re[a^ =n ] (second 
panel). (J,g,(p) = (6,0,3.2). The solid line is obtained when this point is approached from 
the strong coupling regime, and the dashed line is obtained when this point is approached 
from the weak coupling regime. In the second panel, the irregular tails of the dashed line 
are caused by numerical difficulties in the weak coupling regime. 

carries all the crystal momentum. More generally, however, the wave vector of the plane 
wave component is not given strictly by q = n, but by q = n, still reflecting a state in which 
the phonon component carries all the crystal momentum. The latter property suggests that 
outer-zone states on the weak coupling side of this transition are not simply large polaron 
states as we might otherwise characterize them, but mixtures of large polarons essentially 
at rest and "unbound" free phonons. 

The changes in polaron structure that occur as J is varied from small through large values 
at (relatively) fixed 0, can be seen in Figures §a and [| corresponding changes in the polaron 
energy band are illustrated in Figure The particular value of <ft used in this illustration 
was chosen to be significant on an absolute scale, but to fall below both the zone-center and 
zone-edge transitions; this corresponds to a large polaron region in the sense that applies to 
nonlocal coupling. 

At (J, g, (j)) = (0, 0, 1.6), the structure of the polaron state is essentially identical to that 
shown in Fig |^a, with slightly larger amplitudes and weaker decays due to the somewhat 
larger value of <p. The energy band at this point (see Figure 0) has a strongly bimodal 
variation characterized by a negative effective mass at the zone center. These qualitative 
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Figure 8: (3% (first panel) and a* (second panel) for (J, g, • 
and (J,g,(f>) = (1,0,1.6) (dashed line). 



(2,0,1.6) (solid line), 



characteristics of the polaron state and energy band are characteristic of the small J regime. 
At (1, 0, 1.6), on the other hand, the exciton density is essentially completely localized on the 
two central sites, essentially vanishing elsewhere. This coincides with a critcal flattening of 
the polaron energy band at the zone center resulting in a loss of the bimodal band structure 
and a divergence of the effective mass through negative values. With further increases in J, 
e.g. to (2,0, 1.6), the exciton amplitudes again spread in space with a structure dominated 
by J, and the energy band is unimodal with a finite, positive effective mass. 
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Allowing the simultaneous action of local and nonlocal coupling in the presence of finite 
transfer integrals should lead to polaron structures and energy bands that blend the qualities 
seen in Figures |3| and |] ( J = 0; g,(p finite) and Figures || 0, |^, and (g — 0; J, <fi finite). 
Examples of such mixed results for general (J,g,(p) are shown in Figures |10| and [IT| . 



Figure [H] shows how the exciton amplitudes and the phonon displacements vary with J 
for fixed, moderate values of both the local and nonlocal coupling coefficients. 

In the absence of the transfer integral (J = 0), the phonon displacement is mostly 
localized on a single site due to the dominance of local exciton-phonon coupling in this 
case. The effect of nonlocal coupling shows both in the nontrivial spread of the exciton 
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Figure 9: Polaron bands for g = 0, (f> = 1.6, and J = (top, dashed), J = 1 (middle, solid), 
and J = 2 (bottom, dotted), calculated from Toyozawa's Ansatz. 

and phonon amplitudes in space (phonon-assisted transport) and in the alternations in sign 
that persist, though weakened by competition with local coupling. Increasing the transfer 
integral J results in increased spreading that further smoothes both the exciton and phonon 
amplitudes. At J = 3, all amplitudes are positive, and the influence of nonlocal exciton- 
phonon coupling is evident only in the asymmetry of the exciton distribution and its lattice 
distortion. 

In the absence of the transfer integral, Toyozawa's Ansatz and the Munn-Silbey ap- 
proach yield very similar polaron energy bands, though the variational approach based on 
Toyozawa's Ansatz is hereby established as quantitatively superior. It is well to ask, however, 
how superior the present method is, and under what circumstances that difference makes a 
practical difference. The first semi-quantitative conclusion in this regard is that the quan- 
titative discrepancies increase with the transfer integral, confirming one of our expectations 
that motivated this work; the nature of this trend can be seen in a comparison of the lower 
curves of Figure |] (both methods for (0, 1, 1)) with the curves of Figure |ll| (both methods 
for (0.5, 1, 1)). Whether the discrepancies illustrated in Figure |ll| are significant depends on 
the purpose to which a band structure calculation is to be put. The typical "gap" between 
the two results is roughly 0.18 in units of the Einstein frequency, 0.09 in units of the rigid 
lattice energy bandwidth, or about 0.75 in units of the resulting polaron energy band width. 
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Figure 10: Variational parameters for the k = state calculated from Toyozawa's Ansatz. 
First panel: the phonon displacement f3% =0 for g = <p = 1, J = (dashed line), J = 1 (solid 
line) and J = 3 (dotted line); Second panel: the exciton amplitude a* =0 for g = <fi = 1, 
J = (dashed line), J = 1 (solid line) and J = 3 (dotted line). 

For Einstein frequencies typical of molecular phonons, a shift of 0.18 is easily resolved by 
optical spectroscopies, so for such purposes one would certainly want to use the present 
method. On the other hand, the absolute position of the energy band is of less importance 
to other purposes; for example, energy transport in the band regime is more sensitive to the 
shape of the energy band than to its position. The gross shapes of the two energy bands 



in Figure |TT| are similar, though there are signficant differences; for example, the "energy 
barrier" separating the two finite-/? minima is more than twice as large when computed by 
the present method, while the "effective mass" (at k — or at the finite-K minima) is less 
than half that computed by the Munn-Silbey method. 

7 Conclusions 

In this paper we have obtained variational estimates of the ground state energy bands for 
the Holstein Hamiltonian incorporating simultaneous local and nonlocal exciton-phonon cou- 
pling. We have examined the interplay between the transfer integral J and local and nonlocal 
exciton-phonon coupling, with an emphasis on nonlocal coupling effects. 

The most obvious effect of nonlocal exciton-phonon coupling is a bimodal distortion of 
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Figure 11: The polaron bands for (J, g, (f>) = (0.5, 1, 1) at T — 0, calculated from Toyozawa's 
Ansatz (solid line) and the Munn-Silbey approach QTJ (points) . 

the polaron energy band that is most obvious in the limit of small exciton transfer integral 
and small local exciton-phonon coupling. Moreover, since nonlocal coupling is a transfer 
mechanism, a non-trivial polaron energy band remains even when both the transfer integral 
and local coupling vanish. 

The polaron structure induced by local coupling is "site centered", while that induced 
by nonlocal coupling is "bond centered". When local and nonlocal coupling act in concert, 
the microscopic forces associated with each superpose with the result that polaron struc- 
ture is neither completely site-centered nor completely bond centered. The nature of these 
interactions is such, however, that nonlocal coupling must be relatively strong before a bond- 
centered component becomes noticeable against the typically more prominent site-centered 
component driven by local coupling. This latter quality is consistent with findings of other 
analyses of solid-state excimers, where strong nonlocal coupling was found to be an essential 
element in producing the dimeric lattice distortions characteristic of excimers as well as other 
excimeric properties. 

Like polarons born of local coupling only, polarons embracing simultaneous local and 
nonlocal coupling experience self-trapping transitions, manifested in our present calculations 
as in many other approximate treatments by discontinuities in the dependence of some 
polaron properties on system parameters. Unlike the traditional perspective that focusses 
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on the change in the k = polaron state as some "self-trapping line" is crossed in the J — g 
plane, our analysis shows that the traditional notion of self-trapping is a limiting one, the self- 
trapping line corresponding to a one-dimensional boundary of a three-dimensional volume 
in the (J, g, 0) space in which a more general, ^-dependent self-trapping occurs. This more 
general self-trapping phenomenon is associated with polaron bands that are large-polaron- 
like in the inner portion of the Brillouin zone and small-polaron-like in the outer portion 
of the Brillouin zone, the traditional self-trapping corresponding to the limit in which the 
large-polaron-like region vanishes. This characterization of self-trapping continues to hold 
in the J — <f> plane, where the local coupling underlying the traditional self-trapping concept 
no longer exists. In addition to generalizing the traditional self-trapping concept, we find a 
second type of transition phenomenon occurring near the Brillouin zone boundary. Rather 
than constituting a small-to-large polaron type of transition, this phenomeon appears to 
reflect the binding or unbinding of a free phonon. 

While conveniently described as "transitions" , each of these phenomena are best under- 
stood as reflecting rapid but smooth changes in polaron structure that are "marked" by 
discontinuties only because of the approximate nature of our method. For example, there is 
no reason to expect that the "critical point" of one of these transitions has any special mean- 
ing other than marking the point at which the changes in polaron structure first occur too 
rapidly to be accurately represented by our particular calculation method, permitting dis- 
continuities to appear. This quality of self-trapping transitions will receive greater attention 
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For small exciton transfer integrals (J << 1), Toyozawa's Ansatz and the Munn-Silbey 
approach yield very similar polaron energy bands; however, these bands grow increasingly 
dissimilar with increasing transfer integral. In all cases studied, we have found the variational 
approach based on Toyozawa's Ansatz to yield polaron energy bands lower than those of the 
Munn-Silbey method at all k, establishing the variational energy bands as the quantitatively 
superior results. 

Without generalization, the methods used in this paper apply to the limit of zero tem- 
peratures. The Munn-Sibley method, on the other hand, was constructed with the aim of 
maintaining a sound perturbation theory at finite temperatures. Considering the quality of 



REFERENCES 



27 



the comparison between the results of our present method and those of the Munn-Silbey 
method at zero temperature where the latter is not at its best, we would speculate that a 
direct, numerical implementation of the Munn-Silbey method would provide a sound ap- 
proach at finite temperatures, provided that exciton transfer integrals are small relative to 
both local and nonlocal coupling. The variational approach we have employed here needs 
to be generalized in order to reach more reliable conclusions about the finite-temperature 
scenario 
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FIGURE LIST 

Figure 1. The phonon displacement factor calculated from the small polaron Ansatz for 
(J,g,4>) = (0, 1, 0.1). The variation in £™ is small in this case because the nonlocal coupling 
strength (p is weak relative to the local coupling strength g. 

Figure 2. Comparison of the polaron bands between the Munn-Silbey approach Jj]] (points) 
and the small polaron Ansatz (solid line) for (J,g,<p) = (0,1,1,) at T = 0. The small 
polaron Ansatz gives a higher ground state energy and a larger bimodal variation in the 
polaron band. 

Figure 3. Variational parameters for the k = state calculated from Toyozawa's Ansatz. 
Second panel: the exciton amplitude a* =0 (dashed line) and the phonon displacement (3^ =0 
(solid line) for (J, g, (/)) = (0, 0, 1); First panel: the exciton amplitude a^ =0 (dashed line) and 
the phonon displacement j3^ =0 (solid line) for (J,g.<p) = (0,0.3,1). 

Figure 4. Comparison between the polaron band calculated from the Munn-Silbey approach 
Hat T = 0, (diamonds, (0, 1, 1), crosses, (0,0.03, 1) and that from Toyozawa's Ansatz (solid 
line, (0, 1,1), dashed line, (0, 0.03, 1). 

Figure 5. Phase diagram on the J — 4> plane for Toyozawa's Ansatz. The two wedges 
correspond to the discontinuity near the Brillouin zone center (upper wedge) and that near 
the Brillouin zone boundary (lower wedge), respectively. 

Figure 6. Two types of convergent solutions for Re[(3^ =0 ] (first panel) and i?e[a^ =0 ] (second 
panel). (J,g,4>) = (6,0,4). The solid line is obtained when this point is approached from 
the strong coupling regime, and the dashed line is obtained when this point is approached 
from the weak coupling regime. 

Figure 7. Two types of convergent solutions for Re[(3^ =n ] (first panel) and Re[a^ =7T ] (second 
panel). (J,g,(p) = (6,0,3.2). The solid line is obtained when this point is approached from 
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the strong coupling regime, and the dashed line is obtained when this point is approached 
from the weak coupling regime. In the second panel, the irregular tails of the dashed line 
are caused by numerical difficulties in the weak coupling regime. 

Figure 8. (3^ =0 (first panel) and a^ =0 (second panel) for (J, g, <fi) = (2, 0, 1.6) (solid line), and 
(J, g, <fi) = (1, 0, 1.6) (dashed line). 

Figure 9. Polaron bands for g — 0, = 1.6, and J = (top, dashed), J = 1 (middle, solid), 
and J = 2 (bottom, dotted), calculated from Toyozawa's Ansatz. 

Figure 10. Variational parameters for the k = state calculated from Toyozawa's Ansatz. 
First panel: the phonon displacement f3^ =0 for g = <ft = I, J = (dashed line), J = 1 (solid 
line) and J = 3 (dotted line); Second panel: the exciton amplitude a* =0 for g — = 1, 
J = (dashed line), J = 1 (solid line) and J = 3 (dotted line). 

Figure 11. The polaron bands for (J, g, </>) = (0.5, 1, 1) at T = 0, calculated from Toyozawa's 
Ansatz (solid line) and the Munn-Silbey approach [|IJ (points) . 



